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C^ ' Abstract 



Interactions of oppositely charged macroions in aqueous solution give 
rise to intriguing aggregation phenomena, resulting in finite-size, long- 
lived clusters, characterized by a quite narrow size distribution. Particu- 
Q ' larly, the adsorption of highly charged linear polyelectrolytes on oppositely 

Q_) , charged colloidal particles is strongly correlated and some short-range or- 

der arises from competing electrostatic interactions between like-charged 
polymer chains (repulsion) and between polymer chains and particle sur- 
face (attraction). The charge inversion observed for the polyelectrolyte- 
decorated primary particles has been recently explained in terms of corre- 
VQ ' lated adsorption. In these systems, long-lived clusters of polyelectrolyte- 

\f^ . decorated particles form in an interval of concentrations around the in- 

fT^ ' version point. However, the mechanisms that drive the aggregation and 

stabilize, at the different polymer/particle ratios, a well defined size of 
the aggregates are not completely understood. Nor is clear the role that 
fyr\ ' the correlated adsorption plays in the aggregation, although the impor- 

("^ , tance of 'patchy interactions' has been stressed as the possible source of 

attractive interactions between colloidal particles. Different models have 
been introduced to explain the formation of the observed finite-size cluster 
K^ , phase. However a central question still remains unanswered, i.e., whether 

}_^ ' the clusters are true equilibrium or metastable aggregates. To elucidate 

Cu I this point, in this work, we have investigated the effect of the tempera- 

ture on the formation of the clusters. We employed liposomes built up by 
DOTAP lipid interacting with a simple anionic polyion, sodium polyacry- 
late, over an extended concentration range below and over the isoelectric 
condition. Our results show that the aggregation process can be described 
by a thermally-activated mechanism. 



1 Introduction 

Charged colloidal particles in the presence of highly charged linear polyelec- 
trolytes self assemble in long-lived, finite-size mesoscopic aggregates. Despite 
the increasing body of experimental and theoretical work JjiStlSjlHElISlEllHllSl 
[TUl [TT| , the mechanisms that drive the aggregation and stabilize the aggregate 
size are not completely understood. 

The strong and relatively long-ranged electrostatic interactions, the screen- 
ing due to the diffuse ionic layers surrounding the particles, and the simultaneous 
presence of macroions different in geometry and in charge density (the linear 
polyelectrolytes and the bulky particles), make these systems inherently com- 
plex. The formation of the clusters is hence the result of a delicate balance 
between attractions and repulsions characterized by different ranges and also 
differing in their nature (mainly electrostatic ^ but also hydrophobic [121 [5] 
interactions). 

The adsorption of linear polyelectrolytes on the surface of the oppositely 
charged particles occurs in a highly correlated manner. The competing in- 
teractions between polymers and surface (attraction) and between the poly- 
mers (repulsion) results in a locally ordered, correlated adsorption of the chains 
|13[ I14[ I15[ I16[ 117] . This non-uniform distribution gives rise to a non-uniform 
distribution of the surface electric charge, with alternate patches where the 
charge of the polymer or that of the bare particle surface are locally in excess. 
This correlated adsorption is the cause of the counterintuitive phenomenon of 
the 'overcharging', or charge inversion, which occurs when more polyelectrolyte 
adsorbs than is needed to neutralize the original charge of the surface, so that 
the net charge of the polymer-decorated surface inverts its sign. 

From the first observations of the destabilizing effect of polyelectrolytes on 
colloidal suspensions, attractive 'bridging' forces have been proposed as the 
cause of aggregation [TE]. However, when a polyelectrolyte solution and a sus- 
pension of oppositely charged colloidal particles are mixed together, the adsorp- 
tion of the polymer on the particles and the aggregation occur on two different 
timescales. Being favored by the huge surface/volume ratio of the colloid, with 
the surface 'finely distributed' within the whole volume of the host phase, the 
adsorption should be completed in a few seconds from the mixing. Conversely, 
the aggregation process, due to the lower diffusivity of the bulkier colloidal par- 
ticles should occur on a longer timescale. To this reasonable picture has been 
recently given a solid experimental basis [H] . 

As a consequence of the different timescales, the 'primary particles' that are 
involved in the aggregation process are the polyelectrolyte-decorated particles. 
This observation is not to be overlooked, since it suggests a useful shift of the 
perspective to guide the investigation of this phenomenology. Rather than hy- 
pothesizing complex scenarios of 'polyelectrolyte-mediated' interactions, where 
the polyelectrolytes, free in the solution among the particles, drive the aggre- 
gation of the colloid by a mix of osmotic, electrostatic screening and bridging 
effects, one has to deal with a much simpler system of (non-uniformly) polymer- 
coated particles. 



There is much experimental evidence [121 [20] as well as theoretical arguments 
[15l[2Tl[22] in favor of a configuration of the adsorbed polyelectrolyte chains that, 
although with several loops and possibly dangling ends, lie quite flat on the 
particle surface, particularly if the the charge densities on the surface and/or 
on the chains are high [531. Also these considerations suggest that bridging 
mechanism should be excluded as the cause of aggregation, except perhaps for 
very long polyelectrolytes and very small particles. 

The aggregation of charged colloids induced by oppositely charged polyelec- 
trolytes, and the stability of the resulting clusters, appears mainly governed 
by electrostatic forces [231 [13 [H]. While repulsions derive the electrostatic 
diffuse layer overlap, as described in the classical Derjaguin-Landau-Verwey- 
Overbeek theory (DLVO), attractions, besides the ubiquitous van der Waals 
dispersion forces, include an additional non-DLVO attractive force, originating 
from 'patch-charge' heterogeneities |25l[26l[22- In a broad sense, this interac- 
tion can still be viewed as a sort of bridging, since chain configurations span the 
inter-particle region to gain substantial entropy [2H] . 

On the basis of the classical DLVO theory, one expects a stable suspen- 
sion far away from the isoelectric point (IP), and a destabilization, with the 
formation of larger and larger aggregates that eventually precipitate, close to 
the IP. In fact, associated to the progressive reduction of the net charge of the 
polyelectrolyte-decorated particles, an increase of the cluster size is observed. 
The largest size of the aggregates is observed close to the isoelectric point, while 
any further addition of polyelectrolytes promotes a progressive decrease of their 
size. Eventually, when the particle surface is completely saturated by the ad- 
sorbed polyelectrolyte and the overcharging has reached its plateau value, the 
size of the particles in the suspension equals again the size of the primary par- 
ticles ('reentrant condensation') plus a thin layer of adsorbed polymer [31[H1[S]. 
This reentrant condensation has been observed in a variety of colloidal sys- 
tems, ranging from polyelectrolyte-micelle complexes [2], latex particles [29] . 
dendrimers [3D] , phospholipid vesicles (liposomes) [311 [H [S] , to 'hybrid nioso- 
mal' vesicles [T^]. Although the adsorption patterns are further complicated 
by the presence of short-range interactions specific to the different systems, the 
similarities in the observed behavior for such different systems strongly indicate 
that the overall phenomenology is mainly governed by electrostatic interactions, 
arising from double layer overlap (repulsion) and surface charge non-uniformity 
(attraction) . 

DLVO theory does not account for these long-lived finite-size aggregates. By 
itself, the theory only provides a description of the stabilization/destabilization 
behavior. Once the colloid had been destabilized, the size of the clusters evolve 
along dynamics that depends on both particle diffusivity and sticking character- 
istics (diffusion limited or reaction limited cluster aggregation, for example) but, 
in principle, the aggregation eventually results in a phase separation. The for- 
mation of metastable finite-size clusters is not contemplated by the theory, but 
is the result of some different mechanisms. Usually, a metastable cluster phase 
is considered the consequence of a 'frustrated aggregation' due to the decreasing 
diffusivity that progressively reduces the chance of an encounter for the larger 



clusters. Then the metastable cluster phase is 'kinetically stabilized'. However, 
within this theoretical framework, it is difficult to account for the dependence 
of the aggregate size on their net residual charge. 

A different mechanism of stabilization of a cluster phase has been described 
in the case of colloids characterized by the presence of short-range attractive 
and long-range repulsive interactions [321 |331 [Ml [SB I3S]- Different authors 
have hypothesized that stable clusters can result from the equilibrium of long- 
ranged screened electrostatic repulsions and short-range attractions, where the 
like-charged primary particles within the same clusters repel each other even 
though they are far apart and not in a direct contact. However, due to the 
very small screening length in aqueous systems, when compared to the size of 
lipid vesicles, these model does not seem applicable to polyelectrolyte-liposome 
suspensions. 

We have recently shown that the formation of a cluster phase and also the 
observed increase, with the aggregate size, of the apparent cluster-cluster re- 
pulsive interactions as measured by the second virial coefficient [TU], can be 
effectively taken into account the aggregation of the polyelectrolyte-decorated 
particles as a thermally activated process. 

Velegol and Twar [57] have recently developed an analytical model for the 
potential of mean force between non-uniformly charged particles. They have 
shown that a non-uniform charge distribution at the surface of colloidal particles 
gives place to an inter-particle potential that, even in the case of same-sign 
charged particles, has an attractive component. The model is based on the 
Derjaguin approximation and on an extension of the Hogg-Healy-Fuerstenau 
[HHF] model [37]. 

The resulting potential depends on the values of the average electrostatic 
surface potential, <&, and on its standard deviation, a. This inter-particle po- 
tential, for some combinations of the values of $ and a, shows a maximum, 
close to the particles surface, representing an energy barrier that particles have 
to pass for sticking together. The height H of this barrier increases with the 
radius of curvature, R, of the surface of the two approaching particles, i.e., for 
deformable particles (as is the case of polyelectrolyte decorated lipid vesicles), 
with the aggregate size, and depends on the surface potential $, going approx- 
imately a,s H (X R^. By adding to this potential a contribution from van der 
Waals attraction, operating at smaller distance, this picture does not result 
substantially modified. 

As a further test of this picture, in this paper, we investigate the effect of the 
temperature on the aggregation behavior of colloidal charged particles induced 
by addition of oppositely charged polyions. We will show that, as is expected for 
a thermally activated process, the size of the aggregates increases when they are 
formed at higher temperatures. Here, we deal with the complexation between 
an anionic linear polyion, sodium polyacrylate [NaPA], and liposomes built up 
with the cationic lipid di-oleoyl-trimethyl-ammonium-propane [DOTAP] carried 
on at different temperatures, from 5 to 80 °C. Although the range of tempera- 
ture variation is intrinsically limited by the liposome stability and by the aque- 
ous nature of the dispersing medium, a significant effect is observed on cluster 



properties, particularly on their size that, close to the isoelectric point, increases 
appreciably when the temperature is increased from 5 to 80 °C. 

2 Experimental 

2.1 Materials 

Positively charged liposomes were prepared by employing the cationic lipid di- 
oleoyl-trimethyl-ammonium- propane [DOTAP], purchased from Avanti Polar 
Lipids (Alabaster, AL) and used without further purification. The negatively 
charged polyelectrolyte sodium polyacrylate, [—CH2CH{C02Na)—]n [NaPAA], 
with nominal molecular weight 60 kD, was purchased from Polysciences Inc. 
(Warrington, PA) as 25 % aqueous solution. All liposomal samples and poly- 
electrolyte solutions were prepared in Milli-Q water, with electrical conductivity 
less than 1-10~^ mho/cm. 

2.2 Preparation of cationic lipid-polyion complexes 

Liposomes were prepared by dissolution of an appropriate amount of DOTAP 
in methanol-chloroform solution (1:1 vol/ vol). After overnight vacuum roto- 
evaporation of the solvent, the dried lipid film was re-hydrated with Milli-Q 
quality water. The re-hydration process was carried out for 1 hour at a temper- 
ature of 40°C, well above the main phase transition temperature of this lipid 
{Tf « 0°C). In order to form small uni-lamellar vesicles, the lipid solution was 
sonicated for 1 hour at a pulsed power mode, until the solution appeared op- 
tically transparent in white light; the solution was then filtered by means of 
a Millipore 0.4 /xm polycarbonate filter. For all the experiments, liposomes 
were prepared at a concentration of 1.5 mg/ml, corresponding to approximately 
1.5-10^'^ particle/ml (average radius 40 nm) and the solutions were stored at 
4°C. For all the preparations, the size distribution of the DOTAP liposomes 
was log-normal with a mean hydrodynamic radius 40 ± 5 nm and with a poly- 
dispersity of the order of 0.2, as expected for a rather homogeneous particle 
suspension. 

Measurements on liposomes-polyelectrolyte complexes were performed at 
various temperatures, from 5 to 80 °C. Lipoplexes were formed immediately 
before each measurement by mixing equal volumes of the liposome suspension 
and polyion solutions at the appropriate concentration, as elsewhere described 
[H [38] . Before mixing, liposomes and polyelectrolyte solutions were thermostat- 
ted at the desired temperature. The mixed sample was placed in a thermostat- 
ted cell for the measurement of both electrophoretic mobility and size and size 
distribution. Electrophoretic mobility measurements were performed 5 minutes 
after the mixing, immediately followed by size determination. The temperature 
was controlled within ±0.1°C. 



2.3 Electrophoretic mobility and Dynamic light scattering 
measurement s 

Both the electrophoretic mobihty and the size and size distribution of the sus- 
pended particles were measured by means of an integrated apparatus, NanoZe- 
taSizer (Malvern Instruments LTD, UK) equipped with a 5 mW He-Ne laser. In 
this way the electrophoretic parameters and the size of the particles can be mea- 
sured on the same sample and almost simultaneously, reducing the experimental 
uncertainties related to sample preparation, thermal gradients and convective 
movement that, in the presence of large aggregates within the suspension, can 
be considerable. 

The electrophoretic mobility measurements were carried out by means of the 
laser Doppler electrophoresis technique. The mobility u was converted into the 
^-potential using the Smoluchowski relation ( — urj/e, where r/ and e are the 
viscosity and the permittivity of the solvent phase, respectively. 

Size and size distribution of liposomes and polyions- liposome aggregates were 
measured by means of dynamic light scattering (DLS) technique, collecting the 
normalized intensity autocorrelation functions at an angle of 173° and analyzing 
the collected data by using the CONTIN algorithm [39J, in order to obtain the 
decay time distribution of the electric field autocorrelation functions. Decay 
times are used to determine the distribution of the diffusion coefficients D of 
the particles, which in turn can be converted in a distribution of apparent hy- 
drodynamic radii Rh using the StokesEinstein relationship Rh — K bT / QTirjD ^ 
where KbT is the thermal energy and -q the solvent viscosity. The values of the 
radii shown here correspond to the average values on several measurements and 
are obtained from intensity weighted distributions. 

2.4 Simulations 

We simulated the aggregation behavior of a system composed by Np ~ 10000 
particles of initial diameter 2R = 80 nm in a cubic box of volume V with 
packing fraction (f) = AnpR^ /i = 0.01, where p = Np/V is the number density. 
We carried out MC simulation using local metropolis algorithm at T = 298 K. 
Particles interact via a short-range potential defined by eq. [3] The parameters 
of the simulations have been detailed elsewhere [10] . We only recall here that 
to incorporate a Brownian dynamics in the MC algorithm, the zth particle is 
selected with a probability proportional to Rq/Ri, where Rq is the initial radius 
and Ri is the radius of the zth aggregate. Each selected aggregate is then 
moved in each direction by a random quantity (uniformly distributed between 
0.2 nm). Since the aggregation process progressively slows down, simulations 
were interrupted when a plateau in the time dependence of the aggregate average 
radius was reached. 



3 Results and discussion 

We have investigated by means of electrophoretic and dynamic light scattering 
measurements the complexation of cationic DOTAP hposomes in the presence of 
anionic NaPA polyelectrolytes, as a function of the temperature in the interval 
from 5 to 80 °C. 

Preliminarily, the thermal behavior of the liposome dispersion in the ab- 
sence of the polyelectrolyte was investigated. The temperature of the liposome 
dispersion was varied by one step, by immersing the samples, at the initial tem- 
perature of T=25 °C, in a thermostatted bath at the desired final temperature 
and allowing to thermalize for 30 minutes before measuring the size. Follow- 
ing this procedure, no significant variation of liposome size is observed due to 
temperature change. It was found that different thermal protocols induce large 
variations of liposome size, as it was observed for other liposome systems j41j . 
For example, by varying the temperature very slowly, in steps of 1 -^ 2 °C, and 
allowing to thermalize for several minutes after each step, we have observed a 
decrease in the vesicles radius as large as 50% in passing from T=5 to T=:80 °C 
(data not shown). A small, but significant, increase of the ^-potential was ob- 
served at increasing temperature, independently of the procedure employed to 
vary the temperature (data not shown). Taking the advantage of these prelimi- 
nary investigations, only the the procedure of 'rapid immersion' was employed. 

The stability of DOTAP liposome maintained at the desired temperature 
within the interval 5-^80 °C was checked by repeated measurements of size and 
(■-potential. Over a period of 24 hours nor the size nor the (-potential show 
appreciable variations. This high stability of DOTAP liposomes was largely 
expected, since in the whole interval of temperatures investigated the double 
layers of this lipid are in the liquid phase (the main transition temperature of 
this lipid is close to °C) and do not undergo any phase transition that could 
facilitate a lipid restructuring [i^ . 

Analogously, also the complexes that form on mixing liposomes with the 
polyelectrolyte remain stable in time over more than 12 hours, exhibiting a 
mono-modal log-normal distribution when analyzed by CONTIN. Only very 
close to the isoelectric point, large clusters appear, that rapidly precipitate, this 
phenomenon being more pronounced at the higher temperatures investigated. 

As usual 38, 43J , the behavior of the (-potential and hydrodynamic radius 
of complexes is shown as a function of a 'stoichiometric' charge ratio ^, defined 
as the ratio between the total number N~ of the negative charges on the polyion 
chains and the total number iV+ of the positive charges on DOTAP molecules 
in the whole suspension 

iV' _ Cm Mwp 
^ N+ MwM Cd ^ 

where Cm and Cd are the (weight) concentrations of the polyion and DOTAP, 
respectively, and Mwm and Mwd are the molecular weights of the repeating 
unit of the polyion and of DOTAP. This definition of iV^/iV+ considers all the 



charges on the hposomes, both on the outer and on the inner leaflet of the 
vesicle. 

Fig. [1] shows the hydrodynamic diameter 2R and the corresponding ^- 
potential as a function of S, for the liposome-polyelectrolyte complexes, at some 
selected temperatures in the range from 5 to 80 °C. At all the temperatures in- 
vestigated, liposome-polyelectrolyte complexes show a 'reentrant condensation' 
accompanied by a (^-potential inversion, i.e. the inversion of the sign of their net 
electric charge. Close to the inversion point (isoelectric point), the complexes 
reach their maximum size. However, at any given polyelectrolyte-liposome ra- 
tio, the size of the aggregates depends on the temperature, the largest clusters 
being observed at the higher temperatures. This result indicates that aggrega- 
tion is favored by higher temperatures, as it should be in a thermally activated 
process. 

We have recently shown that the observed reentrant condensation and the 
accompanying charge inversion can be both explained in terms of the correlated 
adsorption of polyelectrolytes on the oppositely charged liposomes [3H1 0] . The 
resulting attraction between the polyelectrolyte-decorated particles is due to the 
non- uniform distribution of the net electrostatic charge at their surface [lOl |40] . 

The adsorption of the linear polyelectrolytes on the surface of the oppositely 
charged particles occurs in a highly correlated manner. In fact, the competing 
interactions between polymer and surface (attraction) and between polymers 
(repulsion) result in a locally ordered, correlated adsorption of the polymer 

[HI [Ml [151 [Mini- 

The non-uniform distribution of the adsorbed chains gives also rise to a non- 
uniform distribution of the surface electric charge, with alternate patches where 
the charge of the polymer or the one of the particle surface are locally in excess. 

Velegol and Twar [57] have recently developed an analytical model for the 
potential of mean force between non-uniformly charged particles. They have 
shown that a non-uniform charge distribution at the particle surface results in 
an inter-particle potential that, even in the case of particles that bear a net 
charge of the same sign, has an attractive component. The model is based 
on an extension of the Hogg-Healy-Fuerstenau [HHF] model [37] and on the 
Derjaguin approximation. This approximation holds when the characteristic 
radius of the smaller sphere is much larger than the characteristic length scale 
for the interaction [44], i.e. the Debye screening length. 

By using the Derjaguin approximation, the generic force F{h) between the 
surfaces of two spheres of radii i?i and i?2, at a distance h in terms of the 
potential G(h) that would be observed the two surfaces were infinite planes at 
the same distance /i, can be written as [4^ 

F{h) oc -^^G{h) (2) 

This expression clearly shows that, as a general rule and independently of the 
nature of the inter-particle potential, whenever the Derjaguin approximation 
holds, when their radius increases, the force between two spherical particles 



also increases, toward the limiting force observed for two planes facing each 
other. 

According to Velegol and Twar, the mean force pair interaction potential 
between two spherical particles (A and B) bearing a non-uniformly distributed 
electric charge on their surface can be written, in units of the thermal energy 
fcsT, as 
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were H is the distance between the surfaces of the two approaching parti- 
cles, e the permittivity of the medium and k^^ the Debye screening length. 
This potential depends on the values of the average electrostatic surface poten- 
tials, (^j-potcntial, and on their standard deviations cr,;, combining a net charge- 
dependent monopole term (for Ci 7^ 0)i which is repulsive for like-charged par- 
ticles, and an always attractive multipole term (pi ^ 0) arising from the charge 
heterogeneity. These two components lead up to a global maximum and to 
a potential barrier. The height of the potential barrier that two approaching 
particles must overcome in order to stick together and the separation Umax be- 
tween the particle surface at which this maximum occurs can be evaluated from 
eq. [3] For two identical particles {Ra = ^b — R) we obtain 
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respectively. The height of this barrier increases with the radius of curvature, 
R, of the surface of the two approaching particles, i.e., for deformable particles 
(as is the case of polyelectrolyte decorated lipid vesicles [TU]) with the aggregate 
size. 

At each temperature, the increase of the height of the potential barrier favors 
the formation of finite size aggregates. Within this framework, the observed 
increase of the size of the aggregates at the higher temperatures can be easily 
explained. In fact, as the temperature is increased, the fraction of particles that 
have a sufficient thermal energy to overcome the potential barrier increases. 
Since the barrier height itself increases with the size of the aggregates, larger 
and larger aggregates are stabilized at higher temperatures. 

The aggregation behavior of a system of non-uniformly charged particles 
was recently investigated by the present authors by means of Monte Carlo sim- 
ulations |40J in order to derive the limiting size of the aggregates. Simulations 
qualitatively reproduce the observed trends in different colloidal systems. As 
an example, in fig. [21 it is shown the typical time evolution of the normalized 



mean radius of the clusters that are formed by non-uniformly charged particles 
interacting trough the potential described by equation^ Fig. [2]shows the effect 
on this evolution, for two different average surface electrostatic potential C= 15 
mV and C— H rnV, for a fixed value of the non-uniformity parameter a= 15 
mV. After an initial transient regime the radius of the aggregates increases in a 
similar way for the two considered C-potential values. The growth process slows 
down at long times and the aggregate size reaches a long-time limit value that 
depends on the electrical surface parameters C, and a. For the same value of 
surface inhomogeneity a^ the kinetic arrest is reached at higher values of mean 
cluster size when surface potential C is lower (see inset of Fig. [2|). 

In all our simulations the plateau is reached when the height of the potential 
barrier between the particles is « V^KbT ■ This observation provides a strong 
support to the hypothesis that the dynamics slows down and the arrest of the 
growth observed on the timescale sampled by our simulation is controlled by a 
mechanism of 'thermal stabilization', i.e., the increase of the potential barrier 
with the radius stabilizes a 'limiting size' of the aggregates. Fig. [3] shows the 
comparison between the mean cluster size reached at the kinetic arrest obtained 
by MC simulations (filled circles) and the ones calculated according to eq. H] 
within the framework of Velegol and Twar model, with ^max— 10 ^sT, as a 
function of C, and as a function of a (inset, Fig. [3]). The observed agreement 
supports the evidence for a kinetic arrest at a well defined value of repulsive 
interactions between particles. 

Here, we only report on the effect of the temperature on the size of the 
aggregates. Fig. [4] shows the limiting radius of the aggregates, as a function of 
(■-potential, according to the eq. [U at the different temperatures. In the inset, 
the effect of the charge inhomogeneity parameter a on the mean force potential 
is shown as a function of the minimal distance H between the surfaces of the two 
approaching particles. These behaviors justifies the increase of the size of the 
clusters as the temperature is raised for all the values of the surface potential. 

The experimental results shown in Fig. [1] can be compared in a semi- 
quantitative way with the predictions of the model. Expressing the (-potential 
as a function of the polyelectrolyte-liposome charge ratio ^ (Fig. [U bottom 
panel), the predicted average radius of the aggregates can be calculated as a 
function of ^ from eq. [4] 

Charge inhomogeneity is more difhcult to evaluate. Because of the irregular 
shape and of the high deformability of our aggregates, a direct measurement 
of the a values by electrophoretic rotation |45] appears difficult. However, a 
dependence of cr on ^ can be hypothesizes on the basis of simple considerations. 
a is expected to be close to zero at the beginning of the aggregation process, 
where the aggregates show a (-potential very similar to the one of bare primary 
particles. At these values of the charge ratio, only a few polyion chains adsorb 
individually on the liposome surface and these isolated spots produce negligi- 
ble inhomogeneity on the average. The inhomogeneity of the adsorption layer 
increases with the increase of the polyion content, and a is expected to reach 
its maximum value at the isoelectric point, and then it is expected to decrease 
again as the adsorbed polyion layer becomes more and more compact. This 



10 



behavior can be modeled assuming a Gaussian dependence of a on ^, whose 
width is calculated as the half-height width of the condensation region at the 
different temperatures (Fig. [T] upper panels) , while the amplitude is a free pa- 
rameter to be determined by the fitting procedure. Finally, the value of the 
potential barrier height ^max in eq. 2] is somewhat arbitrarily fixed at 10 i^sT. 
This choice appear justified by MC simulations, showing that for this height a 
dynamical arrest of the aggregation occurs (see, for example. Fig. [2l [3] and ref. 
[10]). With all these assumptions, the eq. |3]can be inverted and the expected 
limiting radius calculated in the whole range of the charge ratio ^ investigated. 

Fig. [5] shows a typical result of this analysis. The diameters measured at 
20 °C (symbols) are compared with the -R(Cj ct) curve calculated on the basis of 
the eq. [4] (continuous line) . An analogous reasonable agreement is found for all 
the temperatures investigated (data not shown). 

Interestingly, the amplitude of the gaussian curve that represents the de- 
pendence of a on ^, i.e. the maximum value assumed by a^ does not vary 
significantly with the temperature (the values being comprised within the inter- 
val 50 ± 5 mV), indicating that inhomogeneity of the surface polyion coverage 
at the isoelectric point is largely independent of temperature. 

It must be noted that the final cluster size has an explicit dependence on 
temperature through the expression of the potential barrier height (eq. S]) . It 
could be interesting to determine if temperature changes affect the cluster size 
according to what is predicted by the explicit dependence on the temperature T, 
or if other effects have to be considered. To this aim, we have plotted together 
all the data at the different temperatures in Fig. [6l shifting the position of 
the aggregation peaks respect to the one at 80 °C. Fig. [6] shows that all the 
experimental data, in the whole range 5-;-80 °C, collapse on a master curve 
when re-scaled, suggesting that the aggregation kinetics of polyion-decorated 
liposome in the temperature range investigated is the same, independently of 
temperature, and it is typical of a thermal activated process. 

This result is expected within the framework of Velegol and Twar model 
[IS]. Considering the dependence of R on (^-potential and a, (eq. S]), once that 
the dependence C(^) and (t(^) are evaluated at each given temperature with the 
aid of experimental data, only the linear dependence on temperature, in the 
potential barrier height, remains. 

At the isoelectric point, where large floes are observed, higher temperatures 
seem to promote a more rapid gross aggregation. The universality of the ag- 
gregation appears to be valid within the condensation region but not at the 
isoelectric point. In fact, data of clusters with almost zero C-potential fail to 
collapse on the master curve, due to large fluctuations of the measured size at 
the isoelectric point, where C-potential has vanishing values. 

Notably, at the lower temperature investigated, 5 °C, on the same exper- 
imental time scale employed in the other measurements, flocculation at the 
isoelectric point is not observed and almost-neutral clusters remain stable. 

The charge ratio ^ where the maximum size of the aggregates is observed 
does not vary significantly at the different temperatures (upper panels, fig. [1]), 
as well as the position of ^-potential sign reversal, corresponding to charge 
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inversion of clusters (bottom panels, Fig. [T|). 

The increase of temperatures (bottom panels, Fig. [1]) promotes a slight 
increase of the ^-potential of the complexes at charge ratios below the isoelectric 
point. This effect is better observed in Fig. [71 where the data are properly 
shifted to make the different positions of the isoelectric point coincident. In 
the regions below and above the isoelectric point, C-potential increases with 
temperature. This small increase can be explained as the result of two combined 
effects. 

First, the effect of counterion condensation on both polyelectrolytes and li- 
posomes must be taken into account. Different theoretical models hypothesize 
that the counterions, on approaching to oppositely charged surface, 'condensed' 
around a highly charged linear polyelectrolyte would not necessarily be released 
because of entropic effects [46l [15l [47]. Recently, these hypotheses have also 
been substantiated by experimental evidence [48l |49l [7] . At the higher temper- 
ature, the increasing entropy of counterions reduces the effects of 'condensation' 
[50] so that the effective charge of the bare surface increases. The second effect 
is connected to the conformational entropy of the adsorbed polyions, that with 
the increase of the polyion content, and hence of the surface coverage, becomes 
relevant. In fact, at low temperatures, polyions tends to assume a more flat 
conformation on liposome surface, thus contributing a stronger screening effect, 
which reflects in the lower values of the <r-potential observed. At higher tem- 
peratures, polyions assume a more disordered conformation with more loops. A 
similar mechanism can be invoked to interpret the behavior observed above the 
isoelectric point. At charge ratios close to 1.5 above the isoelectric condition, a 
sensible increase of the ^-potential (in absolute values) with the temperature is 
observed. At T=80 °C above ^ w 1.5, C-potential values begin to stabilize, as 
if a maximum polyion adsorption were reached. Again, this behavior could be 
connected with the temperature-depending adsorption properties of polyions, 
when, with the increase of the thermal energy, polyion chains are less strongly 
bound to the oppositely charged surface, until a desorption limit is reached [5T] . 

4 Conclusion 

We have investigated the effect of temperature on the reentrant condensation of 
polyelectrolyte-liposome complexes, in a temperature range from 5 to 80 °C. At 
all the temperatures investigated, the observed phenomenology is similar, with 
the formation of polyion-liposome complexes which aggregate in large clusters 
close to the isoelectric point. Temperature does not modify the overall aspect 
of the phenomena. However, significant changes of size and ^-potential values 
of the polyion-decorated liposome aggregates are observed, connected to the 
increase of temperature. 

The observed effects can be interpreted in terms of a thermally activated 
process, where the potential barrier that controls the aggregation arises from 
the combined effect of the electrostatic repulsions, due to the residual net charge 
on the primary particles (the polyelectrolyte-decorated liposomes), and of an 
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attractive term due to the non-uniformity of the adsorbed chains and hence of 
the surface electric potential. 

The observed behavior is interpreted within the framework of Velegol and 
Twar model [57], that takes into account the effect of charge anisotropy on 
aggregating particles. A good agreement is found between calculated and ex- 
perimental values of cluster size, for all the temperatures investigated. More 
interestingly, the analysis according to the Velegol and Twar model evidences 
that, once the dependence of radius on temperature is taken into account, inter- 
actions of oppositely charged polyelectrolytes and liposomes show universality 
properties. In fact, cluster sizes within the reentrant condensation region col- 
lapse on a single master curve, thus excluding the presence of non-electrostatic 
interactions or effects connected with cluster growth rate on the observed ag- 
gregation. 
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Fig. [T]- Reentrant condensation (aggregate diameter 2R, upper panels) and 
charge inversion (^-potential, bottom panels ) of cationic DOTAP liposomes in 
the presence of anionic NaPAA polyelectrolyte. Data are shown as a function 
of the polyion-lipid charge ratio parameter ^ for the different temperatures in- 
vestigated. Lines connecting experimental data point guide eyes only. It has to 
be noted that, at the isoelectrical point, due to the intrinsic instability of the 
system, the data represent the initial size or ^-potential values of flocculating 
aggregates. 

Fig. [2]- Some typical MC-step evolution of normalized mean cluster radius 
< R > / Ro -1. Simulations have been carried out for different values of ^ with a 
constant value of the standard deviation cr=15 mV, o: C=15 mV, o: C,=ll mV. 
The inset shows a magnification of the evolution at longer times of the mean 
size , in the region where kinetical arrest is reached and mean cluster size is 
stable. 

Fig. [3] Comparison between mean cluster size at the kinetic arrest calculated 
according to the eq. H] (continuous hue) with (l>maa. = 10 KbT and obtained by 
MC simulations (filled circles), as a function of the electrostatic parameter C, with 
fixed (7=15 mV. Inset shows mean cluster size according to eq. |4] (continuous 
line) and MC simulations (filled circles), at fixed C=15 mV, as a function of a. 
Each value of cluster size determined by MC simulations is calculated as the 
mean value of the last five points obtained in the MC simulations. 

Fig. [4] - Mean equilibrium radius of clusters at the different investigated 
temperatures as a function of the ^-potential, obtained by means of MC simu- 
lations according to eq. IH with ^^ax= 10 KbT and a= 15 mV (solid hue: 5 
°C, dashed line: 20 °C, dotted line: 40 °C, dash-dot line: 60 °C, dash-dot-dot 
line: 80 °C). Inset shows potential curves as a function of the minimal distance 
between particle surfaces H, at fixed ^-potential (15 mV), calculated for increas- 
ing values of a (sohd hue: 35 mV, dashed line: 30 mV, dotted line: 25 mV, 
dash-dot hne: 20 mV, dash-dot-dot line: 15 mV). 

Fig. [5]- Average hydrodynamic diameter of lipoplex clusters at 20 °C as a 
function of the polyion-lipid charge ratio ^. Line is calculated on the basis of eq. 
[4] within the Velegol and Twar model for inhomogeneously charged interacting 
colloidal particles. 

Fig. [S] - Average hydrodynamic diameter of lipoplex clusters normalized 
with respect to the temperature at which measurement were performed, for the 
various temperatures investigated (D: 5 °C, o: 20 °C, A: 40 °C, y: 60 °C, o: 80 
°C). Cluster size is rescaled with respect to the liposome size and peak position 
is shifted with respect to the position at 80 °C. Errors have been evaluated of 
the order of 10 % of data. Line is to guide eyes only. 
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Fig. [7]- ^-potential curves as a function of the charge ratio £,, for the different 
investigated temperatures (D: 5 °C, o: 20 °C, A: 40 °C, V: 60 °C, o: 80 °C). 
Lines represent hnear fit on data, in regions below and above the isoelectric 
condition. 
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